library(tidyverse)
df<-data.frame(
"period" =factor(rep(c(0, 1), 2), labels=c("before", "after")),
"group" = factor(rep(c(0, 1), each=2), labels=c("control", "treatment")),
"deaths" = c(1349, 1466, 847, 193)
)
model<-lm(deaths ~ period * group , data=df)
tidy(model)|>
select(term, estimate)| term | estimate |
|---|---|
| (Intercept) | 1.35e+03 |
| periodafter | 117 |
| grouptreatment | -502 |
| periodafter:grouptreatment | -771 |
